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Abstract 

In continuation of our perusal of the studies on transverse momentum spectra for the main vari- 
eties of secondaries from a consistent and comprehensive phenomenological approach, wepropose 
to take up here — after a successful completion of reporting in detail the results (Ref. J] in the 



text) on our analyses of the py-spectra of pions — the studies specially on production of kaons, 

> . 

Q\ , protons and antiprotons in several proton- induced and nucleus- involved collisions at high energies. 

00 

The measured data on inclusive cross sections of kaons, protons and antiprotons from the vari- 

o : 

ous published sources have here been assorted first. Next, data on the p^-spectra of the specific 

m 

secondaries produced in PP and PP reactions have been scanned and analyzed with the help of 

43 

Q_i| Hagedorn's model(HM). Thereafter a connector, named here the combinational approach(CA), has 



been constructed and used to analyze the data on p^-spectra of some major category of non-pion 
secondaries produced in nucleus-nucleus (A A/ AB) collisions at high energies. And these exercises 



have, finally, led to the modestly successful interpretations of a wide band of data with the revela- 

?— i ' 

tion of some insightful physical aspects of high energy interactions. The limitations of the approach 
have also been precisely pointed out in the end. 
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I. INTRODUCTION 

The nature of transverse momentum spectra is quite interesting for the following 
reasons: Firstly, in the form of invariant cross section in terms of p?, it offers a ready, reliable 
and very basic observable for both measurements, and for theorization to be attempted. 
Secondly, the pr-spectra throw light on the particle production mechanism on the whole; 
they have also impact and indirect reflections on the predicted features of what is normally 
projected as the signatures of the conjectured quark-gluon plasma(QGP). Thirdly and finally, 
the expressions for invariant p^-spectra lead us to define some other physically significant 
observables in the domain of particle physics, like the average multiplicity, average transverse 
momentum etc. 

After pions, kaons of all varieties (with K ± and K°/K°) constitute the second most abun- 
dant species among the particle secondaries which are produced in multiparticle processes of 
hadron-hadron, hadron-nucleus or nucleus-nucleus interactions at high and very high ener- 
gies. Secondly, kaons are the lightest 'strange' particles for which they have a special status. 
Besides, kaons are supposed to b told above, some strong reflections on what are 

viewed as the QGP diagnostics. Lastly, the kaon-pion ratios and their energy- dependences 
play a very significant role in both the particle physics and non-accelerator physics sectors. 
So, analyzing the kaon production characteristics, especially the p^-spectra, in the light of 
any model assumes a high degree of importance and comes as a challenging task as well. 
Furthermore, the production of proton-antiprotons, the first family-set of baryonic category, 
is also important from various physical considerations. The nature of the ratios P/P has 
a strong bearing on astroparticle physics. Besides, the generalized antiparticle to particle 
ratios are of importance in the study of matter-antimatter symmetry in the universe, since 
they provide the relative abundance of antiparticle pr o duct ions jjj. 

Very recently, spurts of data on the various aspects of particle production from the 
CERN-SPS and the RHIC-BNL, studying high energy interaction properties in the collisions 
between PbPb and AuAu respectively, are available. And providing interpretations for 
both the previous measurements and these oncoming blizzards of new data is really a very 
significant work. 

Our objective here is to try to describe and/or explain the nature of the transverse 
momentum (pr)-spectra of kaons, protons and antiprotons in a large variety of proton- 
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induced and nucleus-induced high energy interactions with the help of a combinational 
approach(CA), called also the grand combination of models (GCM). The present work is 
prompted by our previous successes in analyzing the p r -spectra and rapidity-spectra 0, Q| of 
pions alone in a host of nuclear collisions at high and superhigh energies. In fact, taking up 
the task of interpreting the p^-spectra of a set of non-pion secondaries produced in nuclear 
collisions becomes the logical imperative and the natural follow-up in maintaining a sequence 
of our studies and with the model of our choice in a self-complete way In real terms, while 
proceeding with the present work, we would follow the trajectory of the Ref.Q] here more 
closely than the other one, as in the former the approach to the studies on p^-spectra has 
been rationalized in a better way and much more consistent manner. We have consciously 
and carefully avoided here to adopt operationally any of the known popular/ standard models 
and have also tried to escape, as far as possible, the buzzwords related to them. And we 
have singularly been guided by the results of experimental measurements on the nature of 
Pr-spectra alone in a very dispassionate and objective manner. 

This paper is organized here as follows. In the next section (section II) we give the 
outline of the combinational approach which is to be taken up for this study and the sketch 
of the physical perspective which prompted us to proceed in the stated direction. In section 
III we present the essential steps of the methodology of our work, the results of model- 
based calculations and some brief discussion on the results obtained by this combinational 
approach. The last section is reserved for summing up the conclusions and pointing out the 
shortcomings of the approach adopted here. 



II. THE BASIC APPROACH AND THE PHYSICS PERSPECTIVE 

Following the suggestion of Faesslerpj and the work of Peitzmann .7] and also of Schmidt 
and Schukraft[q], we propose here a generalized empirical relationship between the inclusive 
cross-section for production of K ± , P and P in nucleon(N)-nucleon(N) collision and that 
for nucleus(A)-nucleus(A) collision as given below: 

E % {AB QX) ~ E % (PP QX) ' (1) 

where Q stands for K ± , P, P and 4>(y, px) could be expressed in the factorization form, 

<t>(y, pt) = f(y) g(pt)- 
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While investigating a specific nature of dependence of the two variables(y and px), either 
of them is assumed to remain averaged or with definite values. Speaking in clearer terms, if 
and when p r dependence is studied by experimental group, the rapidity factor is integrated 
over certain limits and is absorbed in the normalization factor. So, the formula turns into 

d 3 n d 3 n 

E— (AB - QX) ~ (A.B)°W E— {PP - QX) , (2) 

The main bulk of work, thus, converges to the making of an appropriate choice of form 
for g{pr). And the necessary choices are to be made on the basis of certain premises and 
physical considerations which do not violate the canons of high energy particle interactions. 

Let us now tentatively propose that the expressions for inclusive cross-section of non-pion 
secondaries in proton-proton scattering at high energies in eqn.(2) could also be chosen in 
the same form as that was suggested for pions by Hagedorn[2|: 

E^ (PP - QX) = & ( 1 + ^)"" , (3) 
dp 6 Po 

where C\ is the normalization constant, and p Q , n are interaction-dependent chosen phe- 

nomenological parameters for which the values are to be obtained by the method of fitting. 

Their y^i-dependences are here proposed to be given by the following formulations: 



MV ~ S) = a + v/i ln (v /i) (4) 



and 

b 

n(Vs) = a + 2 (5) 
In (Vs) 

where a, b, d and b are four constants. The -y/i-dependence of p$ and n would be shown 
later diagrammatically in the text and their data-base would also be indicated. The nature 
and significance of these parameters could be appreciated from the work of Hagedornj^, 
and those of Bielich et al[.l£| and Albrecht et al|ll|. 

The final working formula for the nucleus-nucleus collisions is now being proposed here 
in the form given below: 

E^- {AB -> QX) w C 2 {A.Bf + apT - E—^- {PP -> QX) , (6) 
dp 6 dp 6 

with g{pr) = (e + &Pt — PPt)i wherein this quadratic form of parametrization is 
suggested here tentatively by us for analyzing data and testing its efficacy, if any and is 
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hereafter called De-Bhattacharyya parametrization(DBP). In the above expression C2 is the 
normalization term which has a dependence either on the rapidity or on the rapidity-density 
of the produced secondary; e, a and are constants for a specific set of projectile and target. 

Earlier experimental works 11. 12 . Il^ showed that g(pr) is less than unity in the Pt- 
domain, p? < 1.5 GeV/c. Besides, it was also observed that the parameter e, which gives 
the value of g(pr) at px = 0, is also less than one and this value differs from collision to 
collision. The other two parameters a and essentially determine the nature of curvature 
of g{pr)- However, in the present context, precise determination of e is not possible for the 
following understated reasons: 

(i) To make our point let us recast the expression for (6) in the form given below: 

E^(AB -> QX) « C 2 (A.B) 6 (A.B)^-^) (1 + EL )-« (7) 

Quite obviously, we have adopted here the method of fitting. Now, in eqn.(7) one finds that 
there are two constant terms C2 and e which are neither the coefficients nor the exponent 
terms of any function of the variable, px- And as e is a constant for a specific collision at a 
specific energy, the product of the two terms C2 and (A.B) e appears as just a new constant. 
And, it will not be possible to obtain fit- values simultaneously for two constants of the above 
types through the method of fitting. 

(ii) From eqn.(2) the nature of g(pr) can easily be determined by calculating the ratio 
of the logarithm of the ratios of nuclear-to-PP collision and the logarithm of the product 
AB. Thus, one can measure e from the intercept of g(pr) along y-axis as soon as one gets 
the values of E^% for both AB collision and PP collision at the same cm. energy. In the 
present study we have tried to consider the various collision systems in as many number 
as possible. To do so, we have to consider the data on normalized versions of E%% for 
many collision systems for which clear P^f-data were not available to us. Furthermore, 
from these normalized versions we can/could not extract the appropriate values of as 
the normalization terms, total inclusive cross- sect ions {(T in ) etc., for these collision systems 
cannot always be readily obtained. Besides, it will also not be possible to get readily the data 
on inclusive spectra for PP collisions at all c.m.energies, like e.g., at \fs = 17.8GeV A (c.m. 
energy of Pb + Pb collision) . 

In order to sidetrack these difficulties and also to build-up an escape-route, we have 
concentrated almost wholly here to the values of a and for various collision systems; and 
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the effects of C2 and e have been compressed to a single constant term C3. Hence, the final 
expression becomes, 

E—^-(AB -> QX) « C 3 (A.B) {apT ~^ (1 + ^ )~ n (8) 

with C 3 = C 2 (A.B) e . 

The exponent term apx — (3p\ obviously represents here [gipr) — e] instead of #(pt) alone. 
Thus, after obtaining fit- values of a and /3, if [g(pr) — e] are plotted for various collision 
systems, all the curves would originate from a single point, i.e. origin; and the systems and 
processes are then really comparable. In other words, in this convenient way we could study 
and check the scaling characteristics of g(pr) with respect to the collision systems. 

The expression(8) given above is the physical embodiment of what we have termed to be 
the grand combination of models(GCM) or the combinational approach(CA) that has been 
applied here. Firstly, the results of PP scattering are obtained in the above on the basis 
of an assumption of validity of eqn.(3) for even the non-pion secondaries, though originally 
the Hagedorn's model(HM) was proposed only for pion secondaries. Secondly, the route for 
converting the results for the non-pion secondaries produced in NN to NA or AB collisions is 
built up by the Peitzmann's approach(PA) represented by expression (2). Thirdly, the further 
input is the De-Bhattacharyya parametrization for the nature of the exponent. Thus, the 
GCM or the CA provides the combination of HM, PA and the DBP, and thus constitutes a 
single whole. And this whole is both useful and economical; it is economical in the sense that 
it can accommodate vast wealth of data with only three arbitrary parameters. Compared 
to the theoretical apparatus provided by the standard model(SM)|14|. this number is quite 
low. 

And the choice of this quadratic form is not altogether a coincidence. In dealing with the 
EMC effect in the lepton- nucleus collisions, one of the authors here(SB)Q), made use of 
a polynomial form of A-dependence with the variable ^(Feynman Scaling variable). This 
gives us a clue to make a similar choice with both p? and y(rj) variable(s) in each case sepa- 
rately. In the recent times, De-Bhattacharyya parametrization is being extensively applied 
to interpret the measured data on the various aspects 4] of the particle- nucleus and nucleus- 
nucleus interactions at high energies. In the recent past Hwa et. al. also made use of this 
sort of relationship in a somewhat different context. The underlying physics implications 
of this parametrization stem mainly from the expression(6) which could be identified as a 
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clear mechanism for switch-over of the results obtained for nucleon-nucleon(PP) collision 
to those for nucleus-nucleus interactions at high energies in a direct and straightforward 
manner. The polynomial exponent of the product term on AB takes care of the totality of 
the nuclear effects. 

The individual model(s), seen in a split manner, is(are) certainly not new; but the combi- 
nation with the proposed two-factor quadratic p^-dependence of the exponent (called DBP 
in the text) for AB or AA interactions at high energies is testably a new proposition of- 
fered by us. Besides, the probing of the role of the Hagedorn's model in understanding the 
production of non-pion secondaries in high energy collisions is another virgin feature in our 
work. The testing of this mechanism suggested by us is with the views of (i) presenting 
an integrated approach toward production of various kind of particle-secondaries in high 
energy collisions and (ii) providing a unified outlook to particle production in all varieties 
of particle-particle, particle-nucleus or nucleus-nucleus interactions at high energies. Quite 
understandably, only by intensive analysis of data on several/all sets of collisions somewhat 
successfully such claims could, at all, be made, and/or later be justified. The trends that 
emerge from this extensive data-analysis with this grand combination of models(GCM) has 
been given in the next section. 

Obviouslyin its approach and method, this work is undoubtedly a continuation of one of 
our previous works 4]. But, even at the risk of being a bit repetitive, we have to proceed 
and carry on. Because, without these special efforts at studying on a case-to-case basis for 
each of the specific variety of the secondaries, we cannot reasonably arrive at any definitive 
conclusion and also cannot confirm the physical basis, if any, of our work. And the two points 
made very precisely in the preceding paragraph represent the strong and crucial aspects 
of our physics motivation behind this proposed formalism. A point is to be made. Our 
intention here, is surely not just piling up of large bulk of data in a purposeless way; rather 
we would like to make the fullest possible utilization of them for checking the proposed 
GCM as extensively as possible with the widest range of available data on diverse set of 
reactions. Because, such rigorous and intensive checkings only could lead to meaningful 
and valid conclusions, if any. And, in our opinion, this has helped us to substantiate, to a 
considerable degree, the functional efficacy and strength of the proposed phenomenological 
model which is called here the combinational approach. In fact, the revelation of the potency 
and strength of this phenomenological approach is the other most striking feature which we 
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would surely like to emphasize here. 

It is obvious that there are two phenomenological parameters in g{pr) in expression(8) 
which need to be physically explained and/or identified. In compliance with this condition we 
offer the following physical explanations for the occurrence of all these factors. The particle- 
nucleus or nucleus-nucleus collisions at high energies subsequently gives rise to an expanding 
blob or fireball with rising temperature. In real and concrete terms this stage indicates the 
growing participation of the already-expanded nuclear blob. As temperature increases at 
this stage, the emission of highly energetic secondaries (which are mostly peripheral nucleons 
or baryons) with increasing transverse momentum is perfectly possible. The coefficient a 
addresses this particularity of the natural event; and this is manifested in the enhancement 
of the nuclear contribution with the rise of the transverse momentum. Thereafter, there is 
a turnabout in the state of reality. After the initial fractions of seconds, the earlier-excited 
nuclear matter starts to cool down and there is a clear natural contraction at this stage 
as the system suffers a gradual fall in temperature. Finally, this leads to what one might 
call 'freeze-out' stage, which results in extensive hadronization, especially in production of 
hadrons with very low transverse momentum. In other words, the production of large-pr 
particles at this stage is lowered to a considerable extent. This fact is represented by the 
damping or attenuation term for the production of high-pr particles. The factor {3 with 
negative values takes care of this state of physical reality. Thus the function denoted by 
q{jPt) symbolizes the totality of the features of the expansion-contraction dynamical scenario 
in the after-collision stage. This interpretation is, at present, only of very suggestive nature. 
However, let us make some further clarifications. 

The physical foundation that has here been attempted to be built up is inspired by ther- 
modynamic pictures, whereas the quantitative calculations are based on a sort of pQCD- 
motivated power- law formula represented by eqn.(6). This seems to be somewhat para- 
doxical, because it would be hard to justify the hypothesis of local thermal equilibrium in 
multihadron systems produced by high energy collisions in terms of successive collision of the 
QCD-partons(like quarks and gluons) excited or created in the course of the overall process. 
Except exclusively for central heavy ion collisions, a typical parton can only undergo very 
few interactions before the final-state hadrons 'freeze out', i.e. escape as free particles or res- 
onances. The fact is the hadronic system, before the freeze-out starts, expands a great deal 
- both longitudinally and transversally - while these very few in— tafce p,aceQ. 
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But the number of parton interactions is just one of the several other relevant factors for 
the formation of local equilibrium. Of equal importance is the parton distribution produced 
early in the collision process. This early distribution is supposed to be a superposition of 
collective flow and highly randomized internal motions in each space cell which helps the 
system to achieve a situation close to the equilibrium leading to the appropriate values of 
collective variables including concerned and/or almost concerned quantities. The parameter 
a in expression(6) is somehow related to the measure of the ratio of the net binary collision 
number to the total permissible number among the constituent partons in the pre-freeze 
out expanding stage identified to be a sort of explosive 'detonation' (l2j stage. This is ap- 
proximated by a superposition of collective flow and thermalized internal motion, which is 
a function of hadronic temperature manifested in the behaviour of the average transverse 
momentum. The post freeze-out hadron production scenario is taken care of by the soft 
interaction which is proportional [3, [l^ to the number of participant nucleons, N part , ac- 
cording to almost any variety of wounded nuclear model. The factor (5, we conjecture, might 
have a relationship with the ratio of the actual participating nucleons to the total number 
of maximum allowable(participating) nucleons. In fact, this sort of physical explanations is 
reproduced here from some of our previous works jlQi] . 



III. PROCEDURAL STEPS, RESULTS AND DISCUSSION 

At the very start we study the py-spectra for K ± , P and P inclusive production in PP 
collisions at several energies and try to fit the expression in the formula(3) given in the 
previous section. The graphs are shown in the several diagrams in Fig.l for kaons, protons 
and antiprotons. The fits for the average kaon yield and antiproton production in PP 
reactions have also been obtained(Fig.2) on the basis of Hagedorn's model(eqn.(3)). The 
same could not be done for proton production in PP reactions due to unavailability of the 
measured experimental data in this particular collisions. The obtained fitted-values of the 
parameters, po and n, for different secondaries produced in PP/PP collisions at different 
energies have been displayed in Table-I - Table-V. The graphs for i/s-dependences of p and 
n, as proposed in eqn.(4) and eqn.(5), have been depicted in Fig.3 and the necessary values 
of the parameters, a, d, b and b have been presented in Table- VI. While plotting graphs we 
have obviously assumed that at very high energies the proton-proton and proton-antiproton 
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collisions could be treated at par and with the acceptance of equivalence between each other. 

Hereafter, in studying the nature of pr-spectra in all nucleon-nucleus and nucleus-nucleus 
collisions, the interaction energy in all cases of nucleus-induced reactions is invariably con- 
verted first into the cm. system values, that is expressed in ^/snn, then the values of p and 
n are picked up from the graphical plots drawn already and shown by Fig.3. While analyzing 
nucleus- dependence with expression(8) and trying with the fit parameters C 3 , a and f3, we 
have inserted these extracted values of p and n for PP cross section term occurring in the 
expression(8) given above. 

However, the plots presented in Fig.3 deserve special mention from phenomenological 
points of view. The po and n values for production of the same secondary at various energies 
(^/s- values) show very slow fall-off with energy in these graphs. The nearness of values 
of po/n signals a march towards a steady state of hadronization. Once the natures are 
established by the theoretical procedures adopted for either PP/PP reactions at high or 
very high energies, we would be compelled, for the sake of consistency, to reduce/minimize 
the degree of arbitrariness in choosing values of p and n for nucleon-nucleus or nucleus- 
nucleus collisions. As soon as the values of ^Jsnn is given for a specific interaction, the values 
of po an d n are to be obtained from the theoretical plots of them shown in Fig.3 (a), Fig.3 (b) 
and Fig. 3(c) for different secondaries and/or for different collision(s). The diagrams given 
in Fig. 4 to Fig. 8 do actually provide the nature of collected data and also the theoretical 
descriptions for diverse reactions indicated by the labels in each diagram. The solid curves 
represent the results arrived at on the basis of the new combination of models(NCM) and 
the use of De-Bhattacharyya parametrization(DBP). And the totality of this combination of 
NCM and DBP is called Grand Combination of Models(GCM). And the tables (displayed 
by Table- VII to Table-IX) provide the necessary parameter values that have entered into 
each calculation of the theoretical solid curves in all the diagrams of Fig.4-Fig.8 for various 
reactions. In some cases, due to lack of availability of clear proton data, we have used the 
net proton(P — P) data with the assumption of the near equivalence and as an approximate 
measure. The diagrams presented in Fig.9(a) to Fig.9(c) are full of physical import and 
convey the message of mass number scaling of nuclear reactions. 

The deviations, whenever and wherever if any, from this observation could be attributed, 
in the main, to the following few points: (i) limits of uncertainties in the measurements; (ii) 
some unavoidable approximations in theoretical calculations; and (iii) availability, in some 



10 



cases, of only too sparse and rare data on which reliability is certainly questionable. 

Obviously, we present here the GCM-based analyses of extensive sets of data on some of 
the observables measured by the various experimental groups in high energy particle and 
nuclear physics. The study has, thus, given rise to some crucial observations with revelations 
of following systematic trends: (a) For the non-pion secondaries produced in nucleus-nucleus 
collisions at high energies the nature of fits is not as good as in PP reaction, especially at low 
momentum region, and especially for E^0- vs. pr plots. The reason(s) for this discrepancy 
would be discussed later in the last section, (b) The fit parameters p and n for production 
of secondary kaon(s) and protons are, by ascription, different in the basic PP reactions. 
And this is consistent with the theoretical expectations. The factor po/n has, we recognize, 
a relationship with the slope parameter of the pr-spectra. And the slope parameters for 
the lighter and heavier particles, one knows, are not the same. But, quite interestingly, 
the nature of y^-dependences of p and n values of all three species, kaons, protons and 
antiprotons, as shown in Fig.3(a)-3(b), are qualitatively quite alike. (c) It is quite evident 
from the Tables VII, VIII, and IX that the values of a and (3 are too close. This, in our 
opinion, reflects the basic fact that the nuclear effects are quite finite and certainly not too 
pronounced even with the heaviest of the nuclei and the highest available energies. 

IV. FINAL COMMENTS AND CONCLUSIONS 

In the end let us summarize our main findings from the present study in the following 

way: 

(i) Quite naturally and predictably, the x 2 / n df values shown in Table- VII- IX present, 
at times, sharply contrasting nature, because of either very sparse available data and 
that too for a narrow region of j>r-values, or large ranges of errors and uncertainties 
in the measurements. 

(ii) Still, on the whole, the approach describes modestly well the large bulk of data on 
kaon-antikaon and proton-antiproton production phenomena in high energy nucleon- 
nucleus and nucleus-nucleus collisions. 

(iii) Related with production of K ± and P/P, the cm. energy-dependences of the param- 
eter values occurring in the basic PP reactions are of the same hyperbolic nature as 
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in the cases of pions. This element of observation on similarity aspect of K ± and Pj P 
with pions is certainly a new finding here made by us, as Hagedorn's early work[9j was 
made exclusively on pions. 

(iv) Besides, the entire energy-dependences, even in purely nuclear collisions, are mani- 
fested by only the basic PP interaction; the nuclear geometry exhibits no separate 
energy behaviour. This is quite significant in the sense that it might be a great hint to 
and confirmation of the enormous importance of PP interaction even in understanding 
the nature of heavy ion collisions. 

(v) The observation of the close values of a in various sets of nuclear collisions with varying 
pair of projectile and target and at different energies is another interesting finding. It 
indicates that enhancement due to nuclear contribution could never assume too large 
values, even with the heaviest of the nuclei, though the parameters are uniformly 
non-zero for each and every collision and the values are limited. This does, in effect, 
imply the very finite degree of the nuclear effects. The constraint would appear more 
severe and acute, if the diminutive role of the (3 factor is reckoned with. This might 
be attributed to the large baryon stopping effect in heavy ion collisions. 

Besides, the small positive a-values in any collission represent physically the nu- 
clear enhancement — though certainly not of 'anomalous' nature — called Cronin 
effect, secondly, the negative (3- values might explain the suppressions of productions 
of hadrons at relatively large values of the transverse momenta. 

(vi) Even with some heaviest nuclei, the observations from Fig. 9(a) to Fig.9(c) on the 
validity of a sort of mass-number(A) scaling for production of both strange and non- 
strange variety of hadrons constitute an interesting and probeworthy point. Physically, 
this might signal the limit to the finite degree of impact parameter dependence of 
nuclear collisions. The near-saturation of the parameter value at and after certain 
magnitude of AB indicates the stringent limit on both the number of binary collisions 
and the number of participating nucleons in the colliding nuclei. 

(vii) The reason for observation of departure of the fits from the data on p^-spectra pro- 
duced in nuclear collisions, especially at low region, could be ascribed, in the 
main, to the following reasons. Firstly, because for the nucleus-involved collisions, the 
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measurements in the experiments on purely soft collisions are very difficult, as hard 
collision effects are almost unavoidably present there. Secondly, the measurements 
of the factor E^f through the array of several detectors have always an intrinsic 
simulation-based standard-model (mainly the HIJING version) component, for which 
the obtained data are not the pure products of experimental measurements alone. But 
the present model can, in no way, accommodate such sorts of superposition effects. 

(viii) The property of 'universality' of hadronic, hadronuclear and nuclear reactions at very 
high energies is also vindicated by the approach adopted here. 

The implications of all these are obvious. With the series of such modestly successful 
ventures, the combinational approach attains a plateau and assumes a potential to claim 
the status of an alternative approach to understand functionally some aspects of the heavy 
ion collisions in general, and the nature of ^/rapidity-spectra in NA or AB collisions in 
particular. 

However, we must not gloss over the gross limitations of the present approach which are 
as follows: (i) It is a fact that we could, so far, offer no measure to estimate quantitatively 
the values of a and f3. The hints are, uptil now, of only very qualitative nature, (ii) Quite 
admittedly, the physical explanations outlined here for a and (3 are also of tentative nature; 
they have, still, not been identified concretely with any of the known physical observables in 
the domain of collision dynamics for multiple production of hadrons and of nuclear geometry 
as well, (iii) The calculations for particle-nucleus or nucleus-nucleus interaction cases are 
absolutely dependent on the apriori and data-based knowledge of the production charac- 
teristics of the same particle species in PP reactions at various energies. Without it, the 
method is ineffective and that is a major handicap, (iv) The approach fails to respond to 
the very basic query on actual mechanisms for particle production processes in any collision, 
(v) And just because of it, the approach is insensitive to the charge-state of the specific 
secondary produced in any high energy interaction and that is certainly a great difficulty. 
So, unless these problems could be remedied to a great extent, we accept that we have no 
reasons for complacency. 
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TABLE I: Fit Values of p and n for PP -> (K + + K")/2 + X at different energies 



V NN * ) 

xxti "1" n rpfprPTi 


J. \clc VdllL 
rr~»l 1 lei on— tjnppi ii r*<5 


n. 


P0{yjV V /L) 


7 i 


/AD/ ''^.vjc V / CJ 


X 1 na J 


23[20j 




12 ±3 


2.7 ±0.9 


17 ±4 


0.16 


1.427 


31[20] 


2/cm 0) 


15 ±2 


4±1 


22 ±4 


0.18 


1.031 


45 [20] 


min. bias 


13 ±3 


2.6 ±0.5 


16 ±2 


0.16 


1.276 


53[20] 




15 ±2 


2.4 ±0.2 


15 ±1 


0.16 


1.722 


63[20l 




15 ±1 


1.6 ±0.2 


11 ± 1 


0.15 


1.831 



TABLE II: Fit Values of p and n for PP -> + tf~)/2 + X at different energies 



[with reference ] 


Relevant 
collision-specifics 


Ci 


Po(GeV/c) 


n 


p /n(GeV/c) 


X 2 /ndf 


300 [21] 




5.7 ±0.5 


1.5 ±0.3 


8.2 ±0.5 


0.18 


2.258 


540 [21] 


2/cm < 2.0, 


4.1 ±0.1 


1.3 ±0.4 


8± 2 


0.16 


2.273 


1000[21j 


min. bias 


0.36 ±0.02 


1.50 ±0.06 


7.0 ±0.2 


0.21 


1.971 


1800[21j 




0.09 ±0.01 


1.55 ±0.05 


8±2 


0.19 


3.526 



TABLE III: Fit Values of po and n for PP —> P + X at different energies 



V~s NN (GeV) 
[with reference] 


Relevant 
collision-specifics 


Ci 


Po(GeV/c) 


n 


p /n(GeV/c) 


X 2 /ndf 


23120] 




11 ±2 


19 ±6 


86 ± 15 


0.22 


0.984 


31[20J 




10 ±2 


16 ±5 


70 ± 10 


0.23 


1.012 


45[20J 


Ucm — 0) 


9± 1 


12 ±3 


50 ±9 


0.24 


1.359 


53[20J 


min. bias 


10 ±3 


8±1 


37 ±6 


0.22 


1.022 


63[20J 




10 ±1 


7±3 


34 ± 7 


0.21 


2.106 
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TABLE IV: Fit Values of po and n for PP —>■ P + X at different energies 



^s NN {GeV) 
[with reference ] 


Relevant 
collision-specifics 


Ci 


p (GeV/c) 


n 


PoMGeV/c) 


X /ndf 


93r9f)l 

ZO 1 ZUJ 




ZIZ .^J: 




f\Q -A- Q 


n 99 


1 ns^ 


31[20] 




5.2 ± .5 


13 ±2 


60 ± 11 


0.22 


2.104 


45 [20] 


Dcm — Oj 


7 ±0.4 


11±4 


53 ±8 


0.21 


1.734 


53[20j 


min. bias 


6.6 ±0.7 


10 ±2 


45 ±9 


0.22 


1.430 


63[20] 




9.3 ±0.4 


9±3 


40 ±5 


0.23 


1.226 



TABLE V: Fit Values of po and n for PP — > P + X at different energies 





Relevant 




Ci 


Po(GeV/c) 


n 


p /n(GeV/c) 


X 2 /ndf 


[with reference] 


collision-specifics 












300T21] 






2.6 ±0.2 


7.5 ±0.9 


29 ±3 


0.26 


1.242 


540T21] 


\y cm \ < 2.0 




1.4 ±0.1 


7± 2 


25 ±6 


0.28 


1.775 


1000[21J 


min. bias 




0.13 ±0.03 


7.3 ±0.8 


24 ±3 


0.30 


1.620 


1800[21] 






0.025 ± 0.003 


7.4 ±0.9 


23 ±2 


0.32 


1.483 


TABLE VI: Values of a, 6, d, 6 


Secondaries 




a 




b 




a 


b 


(K+±K-)/2 




1.6 




103 




3.6 


161 


P 




7 




602 




5 


644 


P 




7 




478 




13 


527 
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TABLE VII: Numerical Values of the parameters: a and 

(3 for Kaon production in different high energy collisions 



Collision 


E(GeV) 


Relevant 


c 3 


a 


P 


X 2 /ndf 


[with reference] 




collision-specifics 




(GeV/c)- 1 


(GeV/c)" 2 




P + D[12] 


400 


y = 0, min. bias 


38 ±6 


0.14 ±0.03 


0.04 ±0.01 


2.832 


P + Bef24] 


450 


2.4 < y < 3.5, central 


1.2 ±0.4 


0.26 ±0.03 


0.04 ±0.015 


0.381 


P + S[22] 


200 


0.6 < y < 2.4, central 


6.4 ±0.6 


0.23 ±0.03 


0.03 ±0.01 


1.974 


P + Au[22] 


200 


0.2 < y < 2.0, central 


6.3 ±0.4 


0.21 ±0.03 


0.032 ±0.005 


0.859 


P + PW24] 


450 


2.4 < y < 3.5, central 


1.1 ±0.1 


0.17 ±0.04 


0.04 ±0.01 


0.552 


+ Au[22j 


200,4 


0.2 < y < 2.0, central 


63 ±9 


0.24 ±0.05 


0.03 ±0.01 


1.164 


S + S[22j 


200,4 


0.6 < y < 2.4, central 


72 ± 10 


0.24 ± 0.02 


0.03 ±0.01 


2.088 


S + Pb[24j 


200,4 


2.3 < y < 2.9, central 


28 ±3 


0.25 ±0.03 


0.027 ±0.002 


2.170 


Au + Au[28j 


8450,4 


\Vcm\ < 0.35, min. bias 


17 ±2 


0.15 ±0.03 


0.041 ± 0.005 


1.005 


Pb + Pb[26j 




2.4 < y < 3.5, central 


66 ±2 


0.20 ±0.02 


0.038 ± 0.004 


2.158 
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TABLE VIII: Numerical Values of the parameters: a and 



(3 for Proton 


piOQUCLlOIl 111 


uiiieieni 


nign 


energy collisions 


I tf\ \ 1 l Ql fYn 
vyUlllolUll 


I -J 1 Vjc V J 


Relevant 


c 3 


a 


R 

P 


X 1 na J 


1" rl T'dJT£Jl"£JTI ( ' LA 

W ILil 1 clci cllOt: 1 




collision-specifics 




(GeY/c)- 1 






p j- nri 9i 

r -f- D\LZ\ 


Ann 

41UU 


y = 0, min. bias 


28 ±2 


0.18 ±0.02 


n D98 -i- n nn9 

U.UZo It u.uuz 


9 ^77 
Z.O It 


-r ~r Dei zoj 


410U 


2.3 < y < 2.9, central 


0.57 ±0.02 


0.22 ± 0.04 


U.U4: It U.Ul 


9 

Z.OoU 


r -f- oIZo] 


zuu 


0.5 < y < 3.0, min. bias 


16 ±2 


0.15 ±0.02 


U.U410 It U.UUo 


1 819 


r -(- i\UI Zo] 


zuu 


0.5 < y < 3.0, min. bias 


47 ±6 


0.14 ±0.03 


U.UQ It U.UZ 


1 49^ 


P _L Ph^^l 
r -f- r QIZD] 




2.3 < y < 2.9, central 


0.8 ±0.2 


0.18 ±0.02 


U.U41 It U.ul 


9 

Z. I Do 


n j. An ro^i 


ZUU^i 


0.5 < y < 3.0, min. bias 


102 ±8 


0.16 ±0.02 


n -i- n ri9 

U.UO it U.UZ 


U.OOO 


+ Au[23l 


200,4 


0.5 < y < 3.0, min. bias 


473 ± 63 


0.17 ±0.01 


04 ± 01 


0.866 


S + S[23j 


200,4 


0.5 < y < 3.0, mm. bias 


164 ± 30 


0.19 ± 0.02 


0.03 ± 0.01 


1.247 


S + Ag[23] 


200,4 


0.5 < y < 3.0, min. bias 


450 ± 70 


0.17 ±0.01 


0.03 ± 0.01 


0.992 


S + Au[23] 


200,4 


3.0 < y < 5.0, min. bias 


139 ± 12 


0.17 ±0.03 


0.03 ± 0.01 


2.063 


Au + Au[28j 


8450,4 


\f]crn\ < 0.35, min. bias 


7±2 


0.17 ±0.02 


0.034 ± 0.001 


0.744 


Pb + Pb[26] 


16(L4 


2.3 < y < 2.9, central 


33 ±2 


0.19 ±0.02 


0.022 ± 0.002 


2.672 


TABLE IX: 


Numerical Values of 


the parameters: 


a and (3 


for Antiproton 


production in 


different 


high 


energy collisions 


Collision 


£(GeV) 


Relevant 


c 3 


a 





X 2 /ndf 


[with reference] 




collision-specifics 




(GeV/c)" 1 


(GeV/c)- 2 




P + D[12] 


400 


Ucm = 0, min. bias 


7.2 ±0.6 


0.15 ±0.04 


0.035 ± 0.003 


1.764 


P + Bef25] 


450 


2.3 < y < 2.9, central 


0.16 ±0.03 


0.20 ± 0.03 


0.03 ±0.01 


1.362 


P + Pb[25] 


450 


2.3 < y < 2.9, central 


0.18 ± .04 


0.19 ±0.03 


0.032 ± 0.004 


1.470 


S±S[25] 


200,4 


2.3 < y < 2.9, central 


1.5 ±0.2 


0.19 ±0.02 


0.03 ±0.01 


2.674 


S±Ag[27j 


200,4 


3 < y < 4, central 


8± 1 


0.18 ±0.02 


0.03 ±0.01 


1.851 


S±Au[27j 


200,4 


3 < y < 4, central 


5.7 ±0.6 


0.25 ± 0.03 


0.024 ± 0.003 


1.620 


S±Pb[25j 


200,4 


2.3 < y < 2.9, central 


1.7 ±0.2 


0.23 ± 0.03 


0.022 ± 0.003 


1.767 


Au + Au[28j 


8450,4 


\Vcm\ < 0.35, min. bias 


0.42 ± 0.01 


0.16 ±0.03 


0.034 ± 0.003 


2.081 


Pb + Pb[26] 


160,4 


2.3 < y < 2.9, central 


15 ± 1 


0.26 ± 0.02 


0.024 ± 0.002 


2.838 
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FIG. 1: Plot of vs. pt for various non-pionic secondaries produced in P + P 

collisions at different cm. energies. The various experimental points are taken from 
Ref. j^J. The solid curves give the theoretical fits on the basis of Hagedorn's model(eqn.3). 
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collisions at y/s = 300, 540, 1000 and 1800 GeV. The various experimental points 



are from Ref.j2lJ. The solid curvilinear lines are drawn on the basis of eqn.3. 
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FIG. 3: Values of p and n as a function of cm. energy y/s. Various data points are 
taken from Table-I to Table- V. The solid curves are drawn on the basis of eqn.4 and eqn.5. 
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FIG. 4: Plot of (— )^ as a function of 
transverse momentum, for produc- 
tion of kaons and net protons in PS and 
PAu collisions at E Lab = 200 GeV. The 
various ex per imental points are taken 
from Ref.j3, Q| The solid curves 
depict the GCM-based results(eqn.8). 
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FIG. 5: Plot of Ej§ vs. p T for var- 
ious secondaries produced in P + D 
collisions at E Lah = 400 GeV. The 
various experimental points are taken 
from Ref. [12j. The solid curves pro- 
vide the fits on the basis of eqn.8. 
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FIG. 6: Nature of Invariant spectra for production of K+ + K ; p and P in P + 
Be and P + Pb collisions at -Ex a & = 450 GeV. Various measured data are ob- 
tained from Ref.j^, 13 The solid curvilinear lines are drawn on the basis of GCM. 
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FIG. 7: Plots of multiplicity spectra as function of pt for produc- 
tion of some non-pionic secondaries in different nucleus-nucleus collisions 
at CERN and SPS energies. Various experimental data are taken from 



Ref. 
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The solid curves depict the GCM-based fits(eqn.8). 
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FIG. 8: The nature of invariant spectra for secondary charged kaons, protons and an- 
tiprotons produced in AuAu interactions at RHIC energy (^/s/v n = 130 GeV). The ex- 
perimental data are from Ref. 28]. The solid curves provide the present model-based fits. 
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FIG. 9: Values of a and j3 for various non-pionic secondaries produced in different 
collisions as functions of the product of mass numbers(A_B) of the interacting nuclei. 
The fitted values of a and (3, enlisted in Table-VII to Table-IX, are taken as the 
data points; and are denoted by empty and filled squares respectively. The dashed 
lines give the average values for a K + +K - = 0.21 ± 0.03, ap = 0.17 ± 0.03, ap = 

2 

0.21 ± 0.04, (5 K++K - = 0.035 ± 0.0005, (5 P = 0.037 ± 0.003 and (5 P = 0.030 ± 0.005. 



25 



